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Abstract 

We present the current status of the improved Cascade-Exciton Model (CEM) code 
CEM2k and of the Los Alamos version of the Quark-Gluon String Model code LAQGSM. 
To describe fission and light-fragment (heavier than He4) production, both CEM2k and 
LAQGSM have been merged with the GEM2 code of Furihata. We present some results on 
proton- and deuteron-induced spallation, fission, and fragmentation reactions predicted 
by these extended versions of CEM2k and LAQGSM. We show that merging CEM2k and 
LAQGSM with GEM2 allows us to describe many fission and fragmentation reactions in 
addition to the spallation reactions which are already relatively well described. Neverthe- 
less, the standard version of GEM2 does not provide a completely satisfactory description 
of complex particle spectra, heavy-fragment emission, and spallation yields, and is not yet 
a reliable tool for applications. We conclude that we may choose to use a model similar 
to the GEM2 approach in our codes, but it must be significantly extended and further 
improved. We observe that it is not sufficient to analyze only A and Z distributions of 
the product yields when evaluating this type of model, as is often done in the literature; 
instead it is important to study all the separate isotopic yields as well as the spectra of 
light particles and fragments. 



1 



Introduction 

During recent years, for a number of applications like Accelerator Transmutation of 
nuclear Wastes (ATW), Accelerator Production of Tritium (APT), Rare Isotope Accel- 
erator (RIA), Proton Radiography (PRAD), and others projects, we have developed at 
the Los Alamos National Laboratory an improved version of the Cascade-Exciton Model 
(CEM), contained in the code CEM2k, to describe nucleon-induced reactions at inci- 
dent energies up to 5 GeV PP and the Los Alamos version of the Quark-Gluon String 
Model, realized in the high-energy code LAQGSM j2j, able to describe both particle- and 
nucleus-induced reactions at energies up to about 1 TeV/nucleon. 

In our original motivation, different versions of the CEM and LAQGSM codes were 
developed to reliably describe the yields of spallation products and spectra of secondary 
particles, without a special emphasis on complex-particle and light-fragment emission or 
on fission fragments in reactions with heavy targets. In fact, the initial versions of the 
CEM2k and LAQGSM codes simulate spallation only and do not calculate the process of 
fission, and do not provide fission fragments and a further possible evaporation of particles 
from them. Thus, in simulating the compound stage of a reaction, when these codes 
encounter a fission, they simply tabulate this event (that permits calculation of fission 
cross sections and fissility) and finish the calculation of this event without a subsequent 
treatment of fission fragments. To be able to describe nuclide production in the fission 
region, these codes have to be extended by incorporating a model of high energy fission 
{e.g., in the transport code MCNPX P], where CEM2k and, initially, its precursor, CEM97 
jl], are used, they are supplemented by Atchison's fission model 

Since many nuclear and astrophysical applications require reliable data also on com- 
plex particles (gas production) and light and/or fission fragment production, we addressed 
these questions by further development of CEM2k and LAQGSM codes. We tried differ- 
ent ways of solving these problems and as a first attempt to describe with our codes both 
emission of intermediate-mass fragments heavier than ^He and production of heavy frag- 
ments from fission, we merged CEM2k and LAQGSM with the Generalized Evaporation 
Model (GEM) code by Furihata [3 |H]. We have benchmarked our codes on all proton- 
nucleus and nucleus- nucleus reactions measured recently at GSI (Darmstadt, Germany) 
and on many other different reactions at lower and higher energies measured earlier at 
other laboratories. The size of the present paper allows us to present here only a few 
results, which we choose to be for the GSI measurements on interaction of ^osp]-, ^g^j^is 
with p |9| and d fUJl targets. Results for other reactions may be found in [HI [121 ■ 

CEM2k and LAQGSM Codes 

A detailed description of the initial version of the CEM may be found in Ref. [TBI , 
therefore we outline here only its basic assumptions. The CEM assumes that reactions 
occur in three stages. The first stage is the IntraNuclear Cascade (INC) in which pri- 
mary particles can be re-scattered and produce secondary particles several times prior to 
absorption by or escape from the nucleus. The excited residual nucleus remaining after 
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the cascade determines the particle-hole configuration that is the starting point for the 
preequilibrium stage of the reaction. The subsequent relaxation of the nuclear excitation 
is treated in terms of an improved Modified Exciton Model (MEM) of preequilibrium 
decay followed by the equilibrium evaporative final stage of the reaction. Generally, all 
three stages contribute to experimentally measured outcomes. 

The improved cascade-exciton model in the code CEM2k differs from the older CEM95 
version (which is available free from the NEA/OECD, Paris) |T3| by incorporating new 
approximations for the elementary cross sections used in the cascade, using more pre- 
cise values for nuclear masses and pairing energies, employing a corrected systematics for 
the level-density parameters, adjusting the cross sections for pion absorption on quasi- 
deuteron pairs inside a nucleus, including the Pauli principle in the preequilibrium calcu- 
lation, and improving the calculation of the fission widths. Implementation of significant 
refinements and improvements in the algorithms of many subroutines led to a decrease of 
the computing time by up to a factor of 6 for heavy nuclei, which is very important when 
performing simulations with transport codes. Essentially, CEM2k has a longer cascade 
stage, less preequilibrium emission, and a longer evaporation stage with a higher excitation 
energy, as compared to its precursors CEM97 |4| and CEM95 [14j. Besides the changes 
to CEM97 and CEM95 mentioned above, we also made a number of other improvements 
and refinements, such as: (i) imposing momentum-energy conservation for each simulated 
event (the Monte Carlo algorithm previously used in CEM provides momentum-energy 
conservation only statistically, on the average, but not exactly for the cascade stage of 
each event), (ii) using real binding energies for nucleons at the cascade stage instead of 
the approximation of a constant separation energy of 7 MeV used in previous versions 
of the CEM, (iii) using reduced masses of particles in the calculation of their emission 
widths instead of using the approximation of no recoil used previously, and (iv) a better 
approximation of the total reaction cross sections. On the whole, this set of improvements 
leads to a much better description of particle spectra and yields of residual nuclei and a 
better agreement with available data for a variety of reactions. Details, examples, and 
further references may be found in Refs. pi IT^. 

The Los Alamos version of the Quark-Gluon String Model (LAQGSM) 2J is the next 
generation of the Quark-Gluon String Model (QGSM) by Amelin et al. (see fH] and refer- 
ences therein) and is intended to describe both particle- and nucleus-induced reactions at 
energies up to about 1 TeV/nucleon. The core of the QGSM is built on a time- dependent 
version of the intranuclear cascade model developed at Dubna, often referred in the lit- 
erature simply as the Dubna intranuclear Cascade Model (DCM) (see ^7] and references 
therein). The DCM models interactions of fast cascade particles ("participants") with nu- 
cleon spectators of both the target and projectile nuclei and includes interactions of two 
participants (cascade particles) as well. It uses experimental cross sections (or those cal- 
culated by the Quark-Gluon String Model for energies above 4.5 GeV/nucleon) for these 
elementary interactions to simulate angular and energy distributions of cascade particles, 
also considering the Pauli exclusion principle. When the cascade stage of a reaction is 
completed, QGSM uses the coalescence model described in [T7] to "create" high-energy d, 
t, '^He, and ^He by final state interactions among emitted cascade nucleons, already out- 
side of the colliding nuclei. After calculating the coalescence stage of a reaction, QGSM 
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moves to the description of the last slow stages of the interaction, namely to preequilib- 
rium decay and evaporation, with a possible competition of fission using the standard 
version of the CEM 13 . But if the residual nuclei have atomic numbers with A < 13, 
QGSM uses the Fermi break-up model to calculate their further disintegration instead 
of using the preequilibrium and evaporation models. LAQGSM differs from QGSM by 
replacing the preequilibrium and evaporation parts of QGSM described according to the 
standard CEM ^3] with the new physics from CEM2k ^Ij and has a number of improve- 
ments and refinements in the cascade and Fermi break-up models (in the current version 
of LAQGSM, we use the Fermi break-up model only for A < 12). A detailed description 
of LAQGSM and further references may be found in [2j. 

We have benchmarked CEM2k and LAQGSM against most available experimental 
data and have compared our results with predictions of other current models used by 
the nuclear community. Figure 1 shows examples of calculated neutron spectra from the 
interaction of protons with 208p|j g_g ]^ g GeV compared with experimental data 
PH] ; while Figure 2 gives examples of neutron spectra measured by Nakamura's group (see 
[TU] and references therein) from 560 MeV/nucleon Ar beams on C, Cu, and Pb targets 
compared with our LAQGSM results and predictions by QMD [201 and HIC |21jj from 
Iwata et al. [12]. We see that our codes describe well neutron spectra both for proton- 
and nucleus-nucleus reactions and agree with the data no worse than other models do. 
Similar results are obtained for other reactions for which we found measured data. 

Merging CEM2k and LAQGSM with GEM2 

The Generalized Evaporation Model (GEM) [71 |H] is an extension by Furihata of the 
Dostrovsky et al. [22] evaporation model as implemented in LAHET j23j to include up 
to 66 types of particles and fragments that can be evaporated from an excited compound 
nucleus plus a modification of the version of Atchison's fission model [3 IS] used in LAHET. 
Many of the parameters were adjusted for a better description of fission reactions when 
using it in conjunction with the extended evaporation model. We merged GEM2 (the last 
update of the GEM code) with CEM2k and LAQGSM as follows: we calculate the cascade 
and preequilibrium stages of a reaction with our CEM2k or LAQGSM, then we describe 
the subsequent evaporation of particles and fragments and fission from the remaining 
excited compound nuclei using GEM2. To understand the role of preequilibrium particle 
emission, we performed calculations of all the reactions we tested both with emission 
of preequilibrium particles and without them, i.e., going directly to GEM2 after the 
intranuclear cascade stage of a reaction described by CEM2k or LAQGSM. 

A very detailed description of the GEM, together with a large amount of results 
obtained for many reactions using GEM coupled either with the Bertini or ISABEL INC 
models in LAHET are published by Furihata [301211; many useful details of GEM2 may 
be found also in our paper [TT|. Therefore, we list briefly below only the main features 
of GEM and forward readers interested in details to Refs. [Tj HH 1211 • Furihata did 
not change in GEM the general algorithms used in LAHET to simulate evaporation and 
fission. The decay widths of evaporated particles and fragments are estimated using the 
classical Weisskopf-Ewing statistical theory [23]. 
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Figure 1. Comparison of measured ^8] double differential cross sections of neutrons 
from 0.8 and 1.5 GeV protons on Pb with CEM2k and LAQGSM calculations. 
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Figure 2. Comparison of measured {T3 double differential cross sections of neutrons from 
560 MeV/nucleon Ar beams on C, Cu and Pb with our LAQGSM results and calculations 
by QMD W and HIC ^T] from Iwata et al. jTH]. 
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The new ingredient in GEM in comparison with LAHET which considers evaporation 
of only 6 particles (n, p, d, t, '^He, and "^He) is that Furihata included the possibility of 
evaporation of up to 66 types of particles and fragments (both in the ground and excited 
states) and incorporated into GEM several sets of parameters used to calculate inverse 
cross sections and Coulomb barriers for each ejectile (we use here only default parameters 
of GEM2). The 66 ejectiles considered by GEM are: n, p, d, t, ^'^'^'^Re, ^-^Li, '^'^-^^Be, 
8,10-133^ 10-16C, 14-200, ^^-^ip, ^^-^^Ne, ^'-^'^a, and ^^-^^Mg. 

The fission model used in GEM is based on Atchison's model [21 E] as implemented 
in LAHET often referred in the literature as the Rutherford Appleton Laboratory 
(RAL) model, which is where Atchison developed it. There are two choices of parameters 
for the fission model: one of them is the original parameter set by Atchison [3 E] as 
implemented in LAHET and the other is a parameter set evaluated by Furihata 
[3 IHI , used here default of GEM2. 

The Atchison fission model is designed to only describe fission of nuclei with Z > 70. 
It assumes that fission competes only with neutron emission, i.e., from the widths Tj of 
n, p, d, t, ^He, and ^He, the RAL code calculates the probability of evaporation of any 
particle. When a charged particle is selected to be evaporated, no fission competition 
is taken into account. When a neutron is selected to be evaporated, the code does not 
actually simulate its evaporation, instead it considers that fission may compete, and 
chooses either fission or evaporation of a neutron according to the fission probability Pf. 
This quantity is treated by the RAL code differently for the elements above and below 
Z = 89. The reasons Atchison split the calculation of the fission probability Pf are: (1) 
there is very little experimental information on fission in the region Z = 85 to 88, (2) the 
marked rise in the fission barrier for nuclei with Z'^/A below about 34 (see Fig. 2 in jH]) 
together with the disappearance of asymmetric mass splitting, indicates that a change in 
the character of the fission process occurs. If experimental information were available, a 
split between regions about Z'^/A ^ 34 would more sensible [H]. 

1) 70 < Zj < 88. For fissioning nuclei with 70 < Zj < 88, GEM uses the original 
Atchison calculation of the neutron emission width r„ and fission width Fj to estimate 
the fission probability as 

P = = ^ (1) 

^ Tf + T„. i + r„/r/ 

Atchison uses [3 EI the Weisskopf and Ewing statistical model with an energy- 
independent pre-exponential factor for the level density (see Eq. (6) in JT]) and Dostro- 
vsky's |22| inverse cross section for neutrons and estimates the neutron width F^ as 

F„ = 0.352(1. 68Jo + 1.93A./^Ji + Af^{0.76Ji - 0.05Jo)), (2) 



where Jq and Ji are functions of the level density parameter a„ and Sn[ = 2y a„(£' — Qn — S) 
as 

^ (s„-l)e^" + l ^ ^ (2s,2 - 6s„ + 6)e^" + s2 - 6 

-'0 — 7, ^1^^ -'1 — • 

2a„ 8a„ 

Note that the RAL model uses a fixed value for the level density parameter a„, namely 

an = (A, - l)/8, (3) 
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and this approximation is kept in GEM when calculating the fission probability according 
to Eq. (1), though it differs from the Gilbert-Cameron-Cook-Ignatyuk (GCCI) parame- 
terization (see Eq. (7) in jTT]) used in GEM to calculate particle evaporation widths. The 
fission width for nuclei with 70 < < 88 is calculated in the RAL model and in GEM 
as 

r, ^ (£i^i)£!i±i. (4) 



where sj = 2^af{E — Bj — 5) and the level density parameter in the fission mode aj is 
fitted by Atchison (6j to describe the measured r//r„ as: 

aj = a„(l.08926 + 0.01098(x - 31.08551)2), (5) 

and X = Z^/A. The fission barriers Bf [MeV] are estimated as 

^2 / Z^\^ 

5/ = g„ + 321.2- 16.7-^ + 0.218f-^j . (6) 

Note that neither the angular momentum nor the excitation energy of the nucleus are 
taken into account in the estimate of the fission barriers. 

2) Zj > 89. For heavy fissioning nuclei with Zj > 89 GEM follows the RAL model 
jni El and does not calculate at all the fission width Tj and does not use Eq. (1) to 
estimate the fission probability Pf. Instead, the following semi-empirical expression ob- 
tained by Atchison [3 E] by approximating the experimental values of r„/r/ published 
by Vandenbosch and Huizenga 1211 used to calculate the fission probability: 



log(rjry) = C{Z,){A, - Ao{Z,)), (7) 

where C{Z) and Ao{Z) are constants dependent on the nuclear charge Z only. The values 
of these constants are those used in the current version of LAHET and are tabulated 
in Table 5 of Ref. ^T] (note that some adjustments of these values have been done since 
Atchison's papers 01^1 were published). 

The selection of the mass of the fission fragments depends on whether the fission is 
symmetric or asymmetric. For a pre-fission nucleus with Zf/Ai < 35, only symmetric 
fission is allowed. For Zf/Ai > 35, both symmetric and asymmetric fission are allowed, 
depending on the excitation energy E of the fissioning nucleus. No new parameters were 
determined for asymmetric fission in GEM. 

For nuclei with Zf/Ai > 35, whether the fission is symmetric or not is determined by 
the asymmetric fission probability Pasy 

4870e-o.36i? 

Pn 



asy 



l + 4870e-o-36^' 



For asymmetric fission, the mass of one of the post-fission fragments Ai is selected from 
a Gaussian distribution of mean Af = 140 and width ctm = 6.5. The mass of the second 
fragment is A2 = Ai — Ai. 
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For symmetric fission, Ai is selected from tlie Gaussian distribution of mean Aj = 
Ail 2 and two options for the width om as described in 13 El El- 

The charge distribution of fission fragments is assumed to be a Gaussian distribution 
of mean and width az- Zf is expressed as 

Z, + Z[- 65.5Ai 

Zf = , where Z, = —nr , and / = 1 or 2. 

2 131 + Af/' 

The original Atchison model uses az = 2.0. An investigation by Furihata jH] suggests that 
az = 0.75 provides a better agreement with data; therefore az = 0.75 is used in GEM2 
and in all our calculations. 

The kinetic energy of fission fragments [MeV] is determined by a Gaussian distribution 
with mean ej and width a^j,. The original parameters in the Atchison model are: 

ej = O.USZ^/A]/^ - 11.4, and a,^. = 0.084e/. 
Furihata's parameters in GEM2, which we also use, are: 




0A3lZf/Ay\ for Zf/A^^ < 900, 
OAOAZf/Ay^ + 24.3, for 900 < Zf/AY'"^ < 1800, 



and 




CiiZf/A^' - 1000) + C2, for Zf/A'/' > 1000, 
C2, for Zf/Ay^ < 1000, 



where Ci = 5.70 x 10 ^ and C2 = 86.5. More details may be found in jH]. 

We note that Atchison also has modified his original version using recent data and 
published j2ZI an improved (and more complicated) parameterization for many quantities 
and distributions in his model, but these modifications j2Zj are not yet included either in 
LAHET or in GEM2. 

We have merged the GEM2 code with our CEM2k and LAQGSM, initially keeping 
all the default options in GEM2. We began by concentrating on an analysis of the recent 
GSI measurements in inverse kinematics as the richest and best data set for testing this 
kind of model. As mentioned above, to understand the role of preequilibrium particle 
emission, we performed calculations of all the reactions we tested both taking into account 
preequilibrium particle emission and ignoring it, i.e., going directly to GEM2 after the 
intranuclear cascade stage of a reaction described by CEM2k or LAQGSM. 

If we merge GEM2 with CEM2k without any modifications, the new code does not 
describe correctly the fission cross section (and the yields of fission fragments) whether 
we take into account preequilibrium emission (see the short-dashed line on Fig. 3) or not 
(see the long-dashed line on Fig. 3). Such results were anticipated, as Atchison fitted the 
parameters of his RAL fission model when it was coupled with the Bertini INC j2Hl which 
differs from our INC. In addition, he did not model preequilibrium emission. Therefore, 
the distributions of fissioning nuclei in A, Z, and excitation energy E* simulated by 
Atchison differ significantly from the distributions we get; as a consequence, all the fission 
characteristics are also different. 
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Figure 3. Comparison of tlie experimental mass and charge distributions of the nuchdes 
produced in the reaction p(l GeV) + Pb (circles show data from Tabs. 3 and 4 of Ref. 
[H] and squares - from Fig. 13 of Ref. ^Ul) with different calculations. The dashed lines 
show results found by merging CEM2k with GEM2 without any modifications when pree- 
quilibrium emission is (thin lines) or is not (thick lines) included. Solid lines show results 
from CEM2k+GEM2 with a modified af. thin lines are for the case with preequilibrium 
emission [ofj^^^ j a^^^ = 1.06) and thick lines show the results without preequilibrium 
emission {a^^^^/af"^^ = 0.947). 



Furihata used GEM2 coupled either with the Bertini INC |2H1 or with the ISABEL 
[21] INC code, which also differs from our INC, and did not include preequilibrium particle 
emission. Therefore the real fissioning nuclei simulated by Furihata differ from the ones in 
our simulations, and the parameters adjusted by Furihata to work the best with her INC 
should not be the best for us. To get a good description of the fission cross section (and 
fission-fragment yields) we need to modify at least one parameter in GEM2, namely to 
adjust the level density parameter to get the correct fission cross section (see Eq. (5)), 
in the case of fissioning nuclei with Z < 88 (pre-actinides), and the parameter C{Z) (see 
Eq. (7)) for fissioning nuclei with Z > 88 (actinides). From the dashed lines on Fig. 5 we 
see that we need to enlarge a/ in our code to get a proper fission cross section when we 
include preequilibrium emission (the excitation energy of our fissioning nuclei and their A 
and Z are smaller than provided by the Bertini or ISABEL INC without preequilibrium), 
and we need to decrease in the case without preequilibrium. By increasing af by 1.06 
compared with the original RAL and GEM2 value {a^^^^ /af^^ = 1.06), we are able to 
reproduce correctly with CEM2k+GEM2 the fission cross section for this reaction when 
we take into account preequilibrium emission (below, we label such results as "with Free"). 
In the case with no preequilibrium emission, a proper fission cross section is obtained for 
^CEM I ^RAL _ 947 label such results as "no Free"). We choose these values for a/ 
for all our further CEM2k+GEM2 calculations of this reaction and do not change any 
other parameters. 

The solid lines in Fig. 3 show results with these values of a/. One can see that the 
"no Free" version provide a good description of both the mass and charge distributions 
and agrees better with the data for these characteristics than the "with Free" version 
(that is not true for isotopic distributions of individual elements, as we show below). The 
"with Free" version reproduces correctly the position of the maximum in both A and Z 
distributions and the yields of fission fragments not too far from these maximums, but the 
calculated distributions are narrower than the experimental ones. This is again because 
both Atchison and Furihata fitted their A and Z distributions using models without 
preequilibrium emission, which provide higher values for the excitation energy. A, and Z 
of fissioning nuclei. This means that to get a good description of A and Z distributions 
for fission fragments using GEM2 in CEM2k "with Free" , we would need to modify the A 
and Z distributions of fission fragments in GEM2, making them wider. This would take 
us beyond the scope of the present work and here we do not vary any more parameters 
than we have already discussed. 

Fig. 4 shows the GSI measurements [S] of the A and Z distributions of the kinetic 
energy of products from the same reaction compared with our CEM2k-|-GEM2 calcula- 
tions both with and without preequilibrium emission. Both versions of our calculations 
are in reasonable agreement with the data. 

Mass and charge distributions of the yields or kinetic energies of the nuclides produced 
show only general trends and are not sensitive enough to the details of a reaction. It is 
much more informative to study the characteristics of individual nuclides and particles 
produced in a reaction. Fig. 5 shows a comparison of the experimental data on production 
yields of thirteen separate isotopes with Z lying from 22 to 82 from the same reaction 
measured at GSI ^ with our calculations using both the "with Free" (upper plot) and 
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Figure 4. Comparison of tlie experimental ^Qj mass and charge distributions of the 
spallation-residue kinetic energies of the nuchdes produced in the reaction p(l GeV) + 
Pb (circles) with our CEM2k+GEM2 calculations: "with Prec" results are shown by solid 
lines, "no Prec" results are shown by dashed lines. 



"no Prec" (middle plot) versions. 

The agreement (or disagreement) of our calculations with these data is different from 
what we have for the integral A or Z distributions in Figs. 3 and 4: We see that for the 
isotopes produced in the spallation region (not too far from the target) and for fission 
fragments in the region with the maximum yield, the version "with Prec" agree much 
better with the data than the version "no Prec". Only for production of isotopes at the 
border between spallation and fission and between fission and fragmentation does the 
version "with Prec" underestimates the data, due to too narrow A and Z distributions 
in the simulation of fission fragments, as we discussed previously. The "no Prec" version 
agrees better with the data in these transition regions but are in worse agreement for 
isotopes both in the spallation region and in the middle of the fission region. We conclude 
that if a model agrees well with some A or Z distributions it does not necessarily mean 
that it also describes well production of separate isotopes. In other words, integral A and 
Z distributions are not sensitive enough to develop and test such models, a practice which 
is often used in the literature. 

The lower plot in Fig. 5 shows results of calculations with a version of CEM2k+GEM2 
with reduced preequilibrium emission. We prefer to discus this version together with 
results by LAQGSM+GEM2, and will return to this plot later. 

It is more difficult for any model to describe correctly the energy dependence for the 
production cross sections of different isotopes, i.e., excitation functions. We calculated 
using both the "with Prec" and "no Prec" versions of CEM2k+GEM2 all the excitation 
functions for the same reaction, p + Pb, for proton energies from 10 MeV to 3 GeV and 
compared our results with all available data from our compilation referred to here as T-16 
Library ("T16 Lib") ^30^. Only several typical examples from our comparison are shown 
below. 

Fig. 6 shows two examples of excitation functions for the production of several iso- 
topes in the spallation (first two columns of plots) and fission (the last two columns of 
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Figure 5. Experimental |Qj mass distributions of the cross sections of thirteen isotopes 
with the charge Z from 22 to 82 compared with our CEM2k+GEM2 calculations. "With 
Prec" results are shown on the upper plot, "no Prec" results are shown in the middle, 
while results with reduced preequilibrium emission according to Eq. (9), in the lower one. 
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Figure 6. Excitation functions for the production of ^°^Bi, ^°°T1, ^^^Au, ^^'^Au, ^^^Ir, 
i^^Eu, i24sb, 95Nb, 86Rb, 82Br, 60Co, and ^^Sc from p+^ospb. Results by CEM2k+GEM2 
"with Prec" are shown by sohd hues and "no Prec" by dashed hues. Experimental data 
(filled circles for 208p]^ targets and opaque circles for ""*Pb) are from our LANL compila- 
tion ("T16 Lib") [30j and are available from the authors upon request. 

plots regions). One can see a not too good but still reasonable agreement of both calcula- 
tions with many data (note that most of the data were measured for ""^^Pb targets, while 
our calculations were done for ^^^Pb). We see that merging CEM2k with GEM2 allows 
us to reasonably describe yields of fission fragments, while in the old standard CEM2k 
we do not have any fission fragments and are unable to describe such reactions at all. We 
see that as shown in Fig. 5 for a single proton energy of 1 GeV, the "with Prec" version 
agrees better with the data in the whole energy region both for spallation products and 
for the production of most of the fission fragments. Only on the border between fission 
and fragmentation regions (^^Sc and ^°Co) does the "no Prec" version agree much better 
with the data than the "with Prec" version; the reason for this we have already discussed. 
Similar results were obtained for excitation functions of many other isotopes in the spal- 
lation and fission regions. On the whole, the version "with Prec" reproduces most of the 
experimental excitation functions better that the version "no Prec" . 

In Fig. 7 we show examples of excitation functions for the production of light frag- 
ments (the first two columns of plots), in the fragmentation region, that are produced 
in CEM2k+GEM2 only via evaporation (the contribution to the yield of these isotopes 
from fission or deep spallation is negligible), and of nucleons and complex particles up 
to a (the last two columns of plots). We see that with the "no Prec" version, GEM2 
reproduces correctly the yields of light fragments ^He, ^Li, ^Be, ^'^N, and ^^F, and not so 
well the excitation functions for heavier fragments like ^^Na. With increasing mass of the 
fragment, the calculations progressively underestimate their yields. Note, that in JTj, we 
got very similar results for excitation functions for the p+Au reaction. The version 
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Figure 7. The same as Fig. 6 but for the production of ^^Na, ^^F, ^^N, ''Be, ^Li, ^He, "^He, 
^He, t, d, p, and n. The complex particles up to ^He are produced in CEM2k+GEM2 via 
evaporation and preequilibrium emission; n and p are also produced during the cascade 
stage; fragments heavier than "^He are produced only via evaporation. Data labeled as 
"Segel82" are from Ref. jSIl- 

"with Prec" strongly underestimates the yields of all these fragments, and this is again 
not surprising, as Furihata developed her model and fitted all parameters without taking 
into account preequihbrium processes. 

Undeniably, the parameters determining the yields of evaporated fragments in GEM 
(inverse cross sections and Coulomb barriers) could be adjusted to get a good agreement 
with the data for the yields of light fragments with the version "with Prec" (see, e. g., 
how Furihata and Nakamura addressed this problem in [21] for their version of a code 
without preequilibrium emission). This is not the aim of our present work and we will 
not do this here. Even if we were to do this, we expect in advance to get similar results 
to those we got for the "no Prec" version: It would be possible to describe correctly the 
yields and spectra of light fragments but not of heavy fragments like ^"^Na and ^^Mg. To 
describe such heavy fragments (not only their yields, but also their spectra) the model 
would need to be improved further, by considering other mechanisms for heavy fragment 
production in addition to the evaporation process taken into account by GEM2. 

Finally, the two last columns of plots in Fig. 7 show excitation functions for emission 
of nucleons and complex particles up to a for this reaction. Note that the data for these 
excitation functions are not so extensive and precise as we have for heavier products: many 
data points were obtained by integration (plus extrapolation) of the spectra of particles 
measured only at several angles and only for a limited range of energy. But even from 
a comparison with these sparce and imprecise data we see that the "with Prec" version 
describes these excitation function better than the "no Prec" version, just as we found 
in jTT] for the p+Au reaction. This is an expected result as the high Coulomb barriers 
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for heavy nuclear targets oppose evaporation of low energy charged particles and the 
main contribution to their yields comes from preequilibrium emission from highly excited 
pre-compound nuclei. 

For completeness sake, we show here also an example of results from a calculation 
with the merged CEM2k+GEM2 code of a reaction on an actinide, p(100 MeV) + ^^^U. 
Generally, to get for actinides a proper fission cross section, we need to adjust in GEM2 
the parameters C{Z) (or, also Ao{Z)) in Eq. (7), as they were fitted by Atchison to work 
the best with Bertini's INC and we have in CEM2k our own INC. As mentioned above, 
for actinides, Eq. (1) is not used in GEM2 and aj is not used in any calculations, therefore 
we do not need to adjust a//a„, for fissioning nuclei with Z > 88. We found that for this 
particular reaction, p(100 MeV) + ^^^U, we get with CEM2k+GEM2 a fission cross section 
in agreement with the data without any adjustments of the parameter C{Z) in GEM2, 
i.e., we can use just the default parameters of GEM2. Nevertheless, our results for other 
reactions show that for higher energies of the incident protons or for other target-nuclei, 
the parameter C{Z) has to be fitted to get a correct fission cross section when GEM2 
is coupled either with CEM2k or with LAQGSM. In addition, we should mention that 
for reactions on actinides at intermediate or high energies, the parameter a^^'^^ / aj"^^ 
should also be fitted along with C{Z). In some simulated events several protons can 
be emitted at the cascade and preequilibrium stages of the reaction, as well as at the 
evaporation stage, before the compound nucleus actually fissions (also complex particles 
can be emitted before fission), and the charge of the fissioning nucleus can have Z < 88, 
even when the initial charge of the target has Z > 88. At the same time, for Z < 88, 
due to charge exchange reactions, the charge of the fissioning nucleus may exceed 88, so 
that we would need to fit as well C{Z)'^^'^^ /C{Z)^^^ . This is a peculiarity of treating the 
fission probability Pf differently for the elements above and below Z = 89 in the Atchison 
model. 

Fig. 8 shows mass distributions of products from p(100 MeV) + ^^^U calculated with 
both versions of CEM2k+GEM2 compared to the available experimental data ^] and 
with results by the phenomenological code CYF of Wahl ^33J (short dashed lines). We 
need to mention that these data are not as good for testing and developing models as are 
the GSI data measured in inverse kinematics for the p + Pb reaction discussed above: 
All the data shown in Fig. 8 were obtained by the 7-spectrometry method. Only some 
of the produced isotopes were measured, and most of the data were measured for the 
cumulative yields. To get the "experimental" A-distribution, we summed for each A the 
available data taking care to not sum the individual cross sections already included in 
some cumulative yields; but the resulting A-distribution is still not complete, as many 
isotopes were not measured. This means that some theoretical values can be above the 
experimental data (where some isotopes were not measured) without necessarily implying 
disagreement between calculations and measurements. 

One can see that both the "with Free" and "no Free" versions of CEM2k+GEM2 
describe equally well the mass distribution of the products from this reaction; therefore it 
is not possible to choose between either of the versions from a comparison with these not 
very informative data. We see that results by the phenomenological code CYF of Wahl 
[5^ also agree well with the A-distribution of measured fission products. It is more 
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Figure 8. Comparison of the experimental P2| mass distribution of the nuchdes produced 
in the reaction p(100 MeV) + '^^^U (circles) with calculations by the CEM2k+GEM2 
model when preequilibrium emission is (solid lines) or is not (long dashed lines) included 
and with results by the phenomenological code CYF of Wahl [33^ (short dashed lines). 
Note that measurements were done on ""*U while calculations were performed for ^^^U. 

useful to compare calculations with the yields of individually measured isotopes. In Fig. 
9, we compare our calculations with all cross sections measured by Titarenko [32j for 
each nuclide separately, where we can compare our results with the data (we do not 
include in our comparison the nuclides measured only either in their isomer or ground 
states, as our model does not provide such information: CEM2k+GEM2 provides only 
yields for the sum of isotope production cross sections both in their ground and excited 
states). We see that on the whole, the "with Prec" version agrees better with most of 
the individually measured cross sections than the "no Prec" version and for many of 
the measured isotopes the disagreement is less than a factor of two, with several of our 
calculated points coinciding with the data within the plotting accuracy. Nevertheless, for 
several isotopes like ^^Ga, '^^'^Th'', ^^^Fa% and ^^SNp, we see some big disagreements. For 
comparison, we also show in Fig. 9 calculations by the phenomenological code YIELDX of 
Silberberg, Tsao, and Barghouty and using its 2000 updated version, YIELDX2k [33], 
and with the phenomenological code CYF by Wahl^^ often used in applications. We see 
that these phenomenological systematics fail to describe the production of many isotopes 
from this reaction, indicating that we cannot rely on phenomenological systematics and 
must develop reliable models to be used in applications. 

Very similar results for the studied proton-induced reactions with the ones shown 
above by CEM2k+GEM2 were obtained also by merging LAQGSM with GEM2. We 
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Figure 9. Detailed comparison between experimental j32] (filled squares) and calculated 
cross sections of all individual and cumulative (labeled with a "c") products measured 
for the reaction p(100 MeV) + ^^su. Our CEM2k+GEM2 "with Prec" results are shown 
by connected opaque squares and the "no Prec" results are shown by connected circles. 
Predictions by the phenomenological code CYF [221 shown with connected diamonds, 
and results calculated by the phenomenological systematics YIELDX are shown by 
connected triangles up, and using its 2000 updated version, YIELDX2k j^Sl, by connected 
triangles down. 
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will not discuss such them here; instead, we will apply the merged LAQGSM+GEM2 
code to analyze reactions on the same ^''^Pb-target (projectile, in the case of the inverse 
kinematics of the GSI measurements [IHI), but induced by deuterons. 

A- and Z-distributions of the nuclides produced in the d(l GeV/nucleon) + Pb reac- 
tion calculated with the LAQGSM+GEM2 code are compared in Fig. 10 with the GSI 
data [in]. As in case of proton-induced reactions, we get a correct fission cross section 
only by adjusting the ratio of level density parameters in the fission and evaporation chan- 
nels a//a„ in Eq.(5), or, equivalently, by adjusting afj^^^ j a^^^ , if we use a^^^ provided 
by Eq. (5). We get a correct fission cross section for oP^ ja^^^ = 1.15 in the "with 
Prec" case (short-dashed lines in Fig. 10) and for a'^^^ / aj^^ = 0.97 in the case of "no 
Prec" (long-dashed lines in Fig. 10) when calculating with LAQGSM-I-GEM2 both p- 
and d- induced reactions on Pb at 1 GeV/nucleon. Note that these values differ slightly 
from the ones obtained above for the CEM2k+GEM2 code (see Fig. 3). This can be 
easy understood, as the INC of LAQGSM differs from the INC of CEM2k, therefore the 
nuclei simulated by LAQGSM+GEM2 and CEM2k+GEM2 to really fission after the cas- 
cade and preequilibrium stages of reaction and after evaporation of several particles from 
compound nuclei before fission have slightly different average values of A, Z, and E*. 

The agreement (or disagreement) of the LAQGSM+GEM2 results with these data 
is similar to what we get for the p-|-Pb reaction: the "no Prec" version provides a good 
description of both mass and charge distributions and agrees better with the data for 
these characteristics than the "with Prec" version. These results, together with the ones 
discussed above for p+Pb reactions, as well as results for many other p+A and A-l-A 
reactions measured at GSI at about 1 GeV/nucleon and analyzed with our codes suggest 
us that we need to take into account preequilibrium processes, but we need less emission of 
preequilibrium particles than provided by out standard CEM2k and LAQGSM codes. We 
address this problem following Veselsky It is assumed that the ratio of the number of 
quasiparticles (excitons) n at each preequilibrium reaction stage to the number of excitons 
in equilibrium configuration Ueg, corresponding to the same excitation energy, to be the 
crucial parameter for determining of the probability of preequilibrium emission Ppre- This 
probability for a given preequilibrium reaction stage is evaluated using the formula 



for n < rieq and equal to zero for n > rieq. The basic assumption leading to Eq. (9) is 
that Ppre depends exclusively on the ration n/rieg as can be deduced from the results of 
Bohning j^Tj where the density of particle-hole states is approximately described using 

a Gaussian centered at rieq. The parameter free parameter and no dependence 

on excitation energy is assumed [36^. Our calculations of several reactions using different 
values of apre show that an overall reasonable agreement with available data can be ob- 
tained using apre = 0.4 or 0.5 (see Fig. 11). We choose the fixed value apre = 0.4 both for 
CEM2k+GEM2 and LAQGSM+GEM2 for all our further calculations. 

Results with reduced preequilibrium emission according to Eq. (9) using apre = 0.4 
are shown in Fig. 10 with solid lines (and in the lower plot of Fig. 5 for p+Pb reactions); 
they agree better with the data than either "with Prec" or "no Prec" standard results. 




pre 



(9) 
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Figure 10. Comparison of tlie experimental mass and charge distributions of the nuchdes 
produced in the reaction d(l GeV/nucIeon) + Pb (circles show data from Tabs. 2 and 
3 and squares - from Fig. 13 of Ref. with different calculations. The short-dashed 

lines show results found by merging LAQGSM with GEM2 when preequilibrium emis- 
sion is included and a'j^^^ /a^^^ = 1.15, the long-dashed lines show the results without 
preequilibrium emission for a^^^^ / a^"^^ = 0.97, and the solid lines show the results with 



reduced preequilibrium emission using apre 
details. 



0.4 and a^^^^af^^ = 0.975; see text for 



Fig. 11 shows a comparison of the experimental data on production yields of ten 
separate isotopes with Z lying from 35 to 82 from the same reaction measured at GSI 
HO] with our LAQGSM+GEM2 calculations using both the "with Prec" and "no Prec" 
(upper plot) versions, and when calculated with the version of reduced preequilibrium 
emission according to Eq. (9) (lower plot). One can see that as in the case of the p-|-Pb 
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reaction (see lower plot in Fig. 5), calculations with reduced preequilibrium emission 
agree much better with experimental data than either standard "with Prec" or "no Prec" 
versions. Similar results were obtained for all other p+A and A+A reactions measured 
at GSI and calculated by our codes. Therefore we choose the version with reduced pree- 
quilibrium emission with a fixed value of cXpre of 0.4 as our "standard" version both for 
CEM2k+GEM2 and LAQGSM+GEM2, and all our further calculations will be done in 
this approach. 

Fig. 12 shows a comparison of the experimental JUl A-distribution of the mean 
kinetic energies of the spallation-residue products from the same reaction with different 
calculations by LAQGSM+GEM2. One can see that for this particular reactions all 
calculations agree well with the data, with a slightly better agreement for the version with 
reduced preequilibrium emission. Similar results were obtained for the Z-distribution of 
products from this reaction. 

Finally, for completeness sake, the last four figures (Figs. 13-16) show a comparison 
of all cross sections measured at GSI for the production of nuclides both from spallation 
and fission reactions from interaction of ^°'^Pb beams with p and d targets compared with 
out "standard" LAQGSM+GEM2 results. On can see a very good overall agreement, 
taking into account that we use here our "standard" version of the code and no further 
fitting of any parameters is done. Similar results were obtained with standard versions 
of CEM2k+GEM2 and LAQGSM+GEM2 for many other p+A and A+A reactions mea- 
sured recently at GSI. Such results will be presented in a separate publication. 

Further Work 

Merging the Generalized Evaporation Model code GEM2 by Furihata [3 |H] with 
our CEM2k and LAQGSM codes allows us to describe reasonably well many fission and 
fragmentation reactions in addition to the spallation reactions already described well by 
CEM2k and LAQGSM. But to do so, ones need to provide first a correct description of 
the fission cross sections by fitting in GEM2 the ration a//a„ for pre-actinides and/or 
the parameter C(Z) in Eq. (7) for actinides. This is not a serious problem and we have 
derived already new approximations for a//a„, and for parameter C(Z) in Eq. (7) that 
provide for CEM2k+GEM2 and LAQGSM+GEM2 correct fission cross section calculation 
both for pre-actinides and actinides for a large range of projectile energies. Such new 
approximations will be presented together with examples of our results in a separate 
publication. 

Some reactions, like production of fission fragments at the borders between fission 
and fragmentation or between fission and emission of heavy fragments like Na and Mg 
are poorly described by CEM2k+GEM2 and LAQGSM+GEM2 in their current versions. 
This disagreement does not discourage us; the results of the present work suggest that 
some of the fission and evaporation parameters of GEM2 can be adjusted to get a much 
better description of all reactions. This lends credibility to such an approach. 
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Figure 11. Experimental JUI mass distributions of the cross sections of ten isotopes with 
the charge Z from 35 to 82 compared with our LAQGSM+GEM2 calculations. "With 
Prec" and "no Prec' results are shown on the upper plot, while results with reduced 
preequilibrium emission according to Eq. (9) using a^rec = 0.4 and 0.5 are shown on the 
lower one. 
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Figure 12. Comparison of tlie experimental [TUJ mass distribution of the spallation- 
residue kinetic energies of the nuchdes produced in the reaction d(l GeV/nucleon) + Pb 
(circles) with the same LAQGSM+GEM2 calculations as shown in Fig. 10 . 

There is one more drawback of this approach to be mentioned: considering evapora- 
tion of up to 66 particles in GEM becomes extremely time consuming when calculating 
reactions with heavy targets at high incident energies. But even this disadvantage may 
be mitigated by the performance of modern computers. We have nevertheless some more 
serious doubts about the current version of GEM2 related to its lack of self-consistency, 
e.g.: 

1) using different, not physically related parameterizations for inverse cross sections 
and Coulomb barriers for different particles and fragments; 

2) using different level density parameters for the same compound nuclei when cal- 
culating evaporation and estimating fission probability from the widths of neutron evap- 
oration and fission; 

3) different, and purely phenomenological treatments of fission for pre-actinide and 
actinide nuclei; 

4) not taking into account at all the angular momentum of compound and fissioning 
nuclei; 

5) rough estimations for the fission barriers and level density parameters, etc. 

This means that an approach like GEM2 can in principle be used to describe fission 
and evaporation of particles and fragments heavier that ^He after the INC and preequi- 
librium parts of CEM2k, LAQGSM, and other models, but it should be considerably 
improved striving first to progressively incorporate better physics, and only after that 
adjusting to selected data. 
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Figure 13. Comparison of all measured 9j cross sections of spallation products from 
the reaction 1 GeV/nucleon ^'^'^Pb on p (black circles) with our LAQGSM+GEM2 results 
(open circles). 
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Figure 14. Comparison of all measured [9j cross sections of fission products from the 
reaction 1 GeV/nucleon ^'^'^Pb on p (black circles) with our LAQGSM+GEM2 results 
(open circles). 
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Figure 16. Comparison of all measured ^U] cross sections of fission products from the 
reaction 1 GeV/nucleon ^"^Pb on d (black circles) with our LAQGSM+GEM2 results 
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The results of the present work and from QJJ show that on the whole, the merged 
CEM2k+GEM2 and LAQGSM+GEM2 codes agree better with most of tested experi- 
mental data when we take into account the preequilibrium emission of particles, than 
when we neglect completely preequilibrium processes. But there is still a not so clear 
question here; we have had some indications for many years that GEM accounts for too 
many preequilibrium particles, at least at energies above the pion-production threshold. 
To solve this problem, as a "zero-step" approximation, in our original GEM2k version 
we neglected in the exciton model the transitions that decrease or do not change the num- 
ber of excitons A„ = —2 and A„ = 0, shortening in this way the preequilibrium stage of 
reactions. We do not like this approach as it is somehow arbitrary, "ad hoc" , even though 
this "never come back" approximation is used in some popular codes like LAHET ^2j| 
and FLUKA [3H] . In the present work we removed this arbitrary condition in GEM2k and 
the "with Prec" version takes into account all the preequilibrium transitions A„ = +2, 0, 
and -2, making the preequilibrium stage of a reaction longer and increasing the number 
of emitted preequilibrium particles. The results of the present work indicate to us once 
again that we need to take into account the preequilibrium stage in reactions, but we need 
less particle emission than we calculate at this stage with the standard version of GEM. In 
the present paper we address this problem reducing the number of emitted preequilibrium 
particles following Veselsky 36 , as described above. 

Besides GEM2, we have investigated the well known code GEMINI by Gharity [^2] 
as an alternative way to describe production of various fragments by merging GEMINI 
with both our GEM2k and LAQGSM, and we tested also the thermodynamical fission 
model by Stepanov with its own parameterizations for mass and charge widths, level- 
density parameters, fission barriers, etc., merging it with both GEM2k and LAQGSM to 
describe fission. In addition, we have started to extend GEM2k and LAQGSM and to 
develop our own fission model, as briefly noted in ^T]. The preliminary results we found 
for spallation, fission, and fragmentation products from several reactions we tested so 
far using these approaches are very promising and we will present our results from these 
studies in several separate papers. 
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